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Abstract 

This paper provides primarily an analytical ad hoc -solution for 3- 
dimensional, incompressible Navier-Stokes equations with a suitable 
external force field. The solution turns out to be smooth and inte- 
grable, as there is a gaussian decay of amplitudes present. 

1 Introduction 

This article presents and ad hoc -solution to full Navier-stokes system in 
incompressible case. We begin in a didactic way by presenting a rather trivial 
solution to Navier-Stokes system, then we proceed with a more complex 
situation. The trivial solution is spatially periodic and smooth, the more 
complex solution turns out to be integrable on the whole space, also it is found 
to be a smooth one. The solution is obtained considering a divergence free 
vector field with gaussian decay of amplitudes, and finally adding a suitable 
force field. The mathematical treatment is rather elementary, mainly vector 
calculus. One can refer e.g. [1]. On the details of Navier-Stokes, one can see 
e.g|2]. 
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2 Navier-Stokes equations 



The incompressible Navier-stokes equations are in vector form 

du -> 
— + ^p^iyAu-u-Vu + f (1) 

V-u = (2) 

where u= {ux{x,y, z,t),Uy{x,y, z,t),Uz{x,y, z,t)) is a vector (velocity) field 
u : R-'^x[0,oo) — ^ / = {f:c{x,y, z,t), fy{x,y, z,t), f;,{x,y, z,t)) is a vector 
(force) field / : x [0, oo) — > and p{x, y, z, t) is a scalar (pressure) field 
p : X [0, oo) — > R. > (viscosity) is a constant. 

In terms of Cartesian coordinates, the equations can be stated equivalently 
as 

dt dy dx"^ dy"^ dz"^ ^ 

du.dp_ , f (K\ 

dt ^dz~ dx' + dy^ + dz' ^ dx dy dz + 

duj. duy Q 
dx dy dz 

3 Introductory preliminaries ;solut ion of Navier- 
Stokes with no external force field 

Consider the functions 

Ux{x, y, z, t) = e"*+^(°^+'^^+^") (7) 

Uy{x, y, z, t) = e»*+^(°^+/^2'+^^) (8) 
u,{x,y,z,t) = e'^*+*("^+/^2'+^^) (9) 
with x,y,zEM.,t>0 and a, a, f3,j ^ R. 

We thus have = Uy = = u, equivalently u = {u{x, y, z, t), u{x, y, z, t),u{x, y, z, t)). 
We shall adopt this notation. 

Let us calculate the partial derivatives. We will have the following func- 
tional equations for Ux 



dux 


dux 


dux 


dx 


-"'ay 


dz 


dUy 


dUy 




dx 


-"'ay 


dz 


duz 


duz 


duz 


dx 


-"'ay 


dz 
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^— = laux , — — = iPux , ^— = i-fUx 

ox oy oz 

d Ux 2 ^ n2 ^ ^x 2 

—— = -a Ux , — — = -/3 , — — - = -7 

ox^ oy^ oz^ 
Similarly for Uy 

dUy dUy dUy 

And for Uz 



aa; oy oz 



z 



S^U-, 2 "^^M^ 2 ^^^^ 2 

— = -aUz, = -p Uz , = -7 Uz 

ax^ ay^ az^ 

Prom these we find that 

^ ^ du^ duy duz . .„ . ^ 

V ■ w = h -TT^ + = ^Q:^^x + ipUy + Z7M2 = (10) 

ox oy oz 

As Ux — Uy — Uz — u, we will have 

iu{a + (3 + 'y) = (11) 

Prom which we conclude 

a + p + (12) 
We shall call it the incompressibility condition. 

Let us substitute the partial derivatives into Navier-Stokes equations. We will 
use the notation u = = Uy = Uz- We will obtain 

— — h — ^ = —i'(a'^u + 0^u + 7^m) — iau^ — i^u^ — i'yu'^ (13) 

du dp ,2 ,02 I 2\ • 2 ■ Q 2 ■ 2 ^^ A\ 

— — h 7— = —y\a u + pu + ^uj— lau — ipu — I'yu (14) 
ot oy 

^U dp / 2 , r)2 , 2 \ -2 -a 2 -2 /ir-\ 

— + ^ = —lyya u + pu + ^u)— lau — ipu — i^u (15) 
This is called The Fundamental Functional Equation. 
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Which equal 

^ + ^ = + + f - ^«'(« + /3 + 7) (17) 

Now, as from the incompressibihty condition a + P + j = 0, we will have 
just 

| + | = --.(«^+/3^ + 7> (20) 

1^1 = + + pi) 

This is the fundamental functional equation. Let us turn to calculate the 
time-derivative. As u{x,y,z,t) — e"*+'("^+/^2'+'>'^) , we will have 

du , , 

= «» (22) 

From which we will have the parametric functional equation 

aw+|^ = -i/(Q;2 + /32 + f)« (23) 

ox 

au+^^-u{a'' + P^ + ^^)u (24) 

au + ^ = -u{a' + + 7^)^ (25) 
oz 



This holds, if 

a = -u{a^ + + 7^) and ^ = |^ = |^=0^ p(a;, t/, z,t) ^ c (26) 

c e M. 

This gives the final solution 

u{x, y, z, t) = e-("'+/3'+7')*+^(°^W7^) (27) 

with the restriction ck + /3 + 7 = 0. Dropping the imaginary component, we 
will have 

u{x, y, z, t) = e-'^("'+^'+^')*cos(Q;a: + Py + 'jz) (28) 
with the restriction a + P + j — and 

p{x,y,z,t) = c, c e M (29) 
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4 Analytical solution for Navier-Stokes equa- 
tions 

4.1 The 'grundfunct ions' 

We shall begin with the arbitrary choice of 'grundf unctions', as we shall call 
them. They are completely ad hoc in their nature. 

u^{x, y, z, t) = a{Biy - y){Ciz - ^)e«*-|(^'+2''+-')+i(^^'+S2''+t^-') (30) 

Uy{x,y,z,t) = p{Atx - x){Ciz - ^)e«*-|(^'+2''+-')+f (^^'+^2^'+^-') (31) 
u,{x, y, z, t) = ^{Aix - x){Biy - y)e^t-\(^'WW)+i{Ax^+By^+Cz^) (33) 
with a, A, B, C,a, & C and i is the imaginary unit. 

For things to come, we calculate straightforwardly 

u^Uy = aP{Aix - x){Biy - y){Ciz - zfe^<^t-i^'W+z^)+i{Ax^+By^+Cz^) (33) 

M^M, = a^{Aix - x){Biy - yf{Ciz - z)e^'^t-{^'W+z^)+i{Ax^+By^+Cz^) (34) 

UyU, - (5-i{Aix - xf{Biy - y){Ciz - z)^'^i-i''^WW)^i{Ax-^^By^^Cz-^) ^35) 

We note immediately that the divergence of such vector field of grund- 
f unctions vanishes, as 

V ■ M = {^Aix — x)ux + {Biy — y)uy + {Ciz — z)uz — 0\/x, y, z, t (36) 

if and only if 

a + /3 + 7 = (37) 
This is the incompressibility condition. 



4.2 Derivatives of grundfunct ions 

Let us calculate the partial derivatives. We will have the following functional 
equations for Ux 

— = {Aix - x)Ux 

, . , Ux 

^— = [Biy - y)ux H 

oy y 
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^— = [Ciz - z)ux H 

az z 

The second partial derivatives 

= {Aix - xfu^ + {Ai - l)u^ 

= {Biy - yfux + 3(Si - l)itj 
= {Ciz — z^Ux + 3(Cz — l)ux 



The cross-terms 



Uy— = [Biy - y)UxUy H 

oy y 

Uz^- = (<->«^ - Z)UxU^ + 



X "'Z 



So sum of cross-terms will be 

Ux{{Aix - x)ux + {Biy - y)uy + ^ + {Ciz - z)uz + — ) = + (38) 

y z y z 

And for Uy 

^^{Mx-x)uy + ^ 

dUy . . 

- = {Biy - y)uy 



dy 

^ = {Ciz - z)Uy + ^ 

oz z 
The second partial derivatives 

Y — {Aix — x^Uy + 3{Ai — l)uy 



dx^ 



^ - {Biy - yfuy + {Bi - l)uy 



dy"^ 



= {Ciz - zfUy + 3(CZ - l)Uy 
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The cross terms 



'^x-g^ = Ux{{Aix - x)uy + 
""y-gf = iBiy - y)ul 

dUy , Uy 



^ dz 

So sum of cross-terms will be 



= Uz{(Ciz — z)uy + — ) 



Uy{{Aix - x)u, + {Biy -y)uy + ^ + {Ciz - z)u, ^^^^^^^ (39) 

X Z X z 



And for u 



z 



-;r— = [AlX - X)Uz H 

OX X 

— = [Biy - y)u^ H 

dy y 



dz 

The second partial derivatives 



{Ciz — z)uz 



— = (Aix - xfuz + 3(Ai - l)uz 



dz^ 



{Biy - yfuz + 3(Si - l)^^ 
: {Ciz — z^Ux + {Ci — l)ux 



The cross-terms 



Ux^w^ = Ux{{Aix - x)uz + — ) 
ax X 

Ou U- 
Uy-Q^ = Uy{{Biy - y)u^ + 

Uz-g^ = {ClZ - Z)U^ 

So sum of cross-terms will be 

u,{{Aix-x)ux + {Biy-y)uy + — + {Ciz-z)u,+ '^) ^ !f£!^ + !f£% (40) 

X y X y 
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4.3 Navier-Stokes equations with the choice of grund- 
functions 

The right side of the Navier-Stokes -system will become 

?/ ?/ 11 II 

i/((Aix-x)2 + (Ai-l) + (Si2/-2/)2+3(S^-l) + («^-0)2-i-3(Cz-l))M^-(^+-^)+/^ 



z/((A*x-x)^+3(Ai-l) + (S^y-|/)^+(Sz-l)+(Cz^-z)^+3(Cz-l))M2,-(^^+^^)+/2 



y 


z 


(41) 




UyUx 1 




X 


z 


(42) 









X y 
(43) 

The time derivatives will be obviously 

r = (44) 



dt 
du 



dt 

du. 



^ = auy (45) 
= aUz (46) 



dt 

Therefore we naturally choose, 

a^v{{A + B + C)i-^) (47) 

As we want a to be a constant. Therefore the grundfunctions will become 

Ux{x, y, z, t) = a{Biy - y){Ciz - ^)e-((^+^+<^)^-3)*-|(^'+2/'+-')+t(^-' W+c.^) 

(48) 

Uy{x, y, Z, t) = (5{Aix - x){Ciz - ;,y{{A+B+C)i-Z)t-\[.'^+y-W)n^A.'^+By^+Cz^) 

(49) 

u,{x, y, z, t) = ^{Aix - x){Biy - yyiiA+B+c)i-3)t-Ux-W+z')Hi^x-+By'+Cz^) 

(50) 

One can see that the the field is rotating with respect to time due to the 
factor {A + B + C)i. 

Therefore the remaining pressure gradient will have components 

^ = K(^^^-^)'+(5^y-2/)'+2(i^^-l)+(^7^^-^)'+2(C^-l))u,-(^+^)+/, 

^ (51)"^ 



X 



^ = z.((Aix-x)2+2(^z-l)+(5zy-y)2+(Cz^-z)2+2(Cz-l))«,-(^+^)+/, 

^ (52)' 

1^ = ^.((Aix-a;)^+2(Ai-l)+(Biy-y)V2(Bi-l)+(«^-^)^)i.,-(^+^)+/, 

X y 

(53) 

We do not know in general, whether such pressure field exists, instead, we 
construct an explicit solution. 

4.4 Explicit solution with a nontrivial external force 
field 

We will show an explicit solution with the aid of an external force field. We 
consider an external force field 

f-if.JyJz) (54) 

We assume it is smooth and intcgrable. For simplicity we will take A = B = 
C = 1, that is, the frequencies arc all of unity. Now consider the pressure 
gradient. 

1^ = u{{ix-xfu^+{iy-yfu^+2{i-l)u^+{iz-zyu^+2{i-l)u^)-{^^^+^^^ 
ox y z 

(55) 

— = i'{{ix-x)'^Uy+2{i-l)uy+{iy-yYuy+{iz-z)'^Uy+2{i-l)uy)-{-^-\ — —)+fy 
oy X z 

(56) 

^ = u{{tx-x)\,+2{t-l)uMw-yfu,+2{t-l)u,+{tz-z)\,)-{^+^)+f, 



oz X y 

(57) 

Let us multiply the first equation by {xi — x) , the second by {yi — y) and 
the third by {zi — z), and sum the three equations with each other. Using 
the property of zero-divergence, we will get just 

-x(^-l)(^+^-/,)-y(^-l)(^+^-/,)-^(^-l)(^+^-/,) = 
y z X z X y 

(58) 

+ + (59) 

Let us choose backwards so that 

d^^ u^^u^ (60) 
ox y z 
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dp 



dy X z 



As we have 



dz X y 

—AaPxz^e' 



y 

z 



-Apayz'^e'' 
—Aa'yzy'^e'' 



X 

X 



z 

With g — e 



u{6i-6)t-(x^ +1/2 +z'^)+i(x^ +1/2 +z2 ) 



^ = Aa{pz^ + -fy^)xg + 
^ = 4/3(q;2;' + 7x')|/^ + fy 

Now suppose we have an external force field /, such that 

= -AaiPz'^ + -fy'^)xg + 4xg 

fy = -4/?(a^^ + 'yx'^)yg + Ayg 
f, = -A^iay^ + I3x^)zg + 4.zg 
Then the pressure gradient will be just 

dp . 
dp . 
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This can be integrated easily to obtain 

p{x, y, z, t) = -(1 + ^)e-(6^-6)t-(^'+2/'+-')+^(^'+2/'+-') + constant (78) 

4.5 Conclusion 

We have provided an analytical solution to the full Navier-Stokes equations 
that is smooth and integrable with the aid of an external force field. Now 
one can proceed to extend this solution to more complex situations. What 
kind of frequency combinations are possible? Are there nontrivial solutions 
without an external force field? 
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